Evidence for Relativistic Outflows in Narrow-line Seyfert 1 Galaxies 
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^ ; ABSTRACT 

o : 

Q"^ . We report the observation of features near 1 keV in the ASCA spectra from three 

fH I "Narrow Line Seyfert 1" (NLSl) galaxies. We interpret these as oxygen absorption in a 

highly relativistic outflow. If interpreted as absorption edges, the implied velocities are 
O ' 0.2-0.3 c, near the limit predicted by "line-locking" radiative acceleration. If instead 

' interpreted as broad absorption lines, the implied velocities are ~ 0.57 c, interestingly 

near the velocity of particles in the last stable orbit around a Kerr black hole, although 
J> ■ a physical interpretation of this is not obvious. The features are reminiscent of the UV 

^ . absorption lines seen in broad absorption line quasars (BALQSOs), but with larger 

^ ! velocities, and we note the remarkable similarities in the optical emission line and 

broad band properties of NLSls and low-ionization BALQSOs. 



Subject headings: galaxies: individual (IRAS 13224-3809, IH 0707-495, PC 1404+226) 
- X-rays: galaxies - galaxies: active 
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1. Introduction 

Narrow-line Scyfert 1 galaxies (NLSls) are defined by their optical line properties (e.g. 
Goodrich 1989): (i.) the Balmer lines are only slightly broader than the forbidden lines 
(H/? FWHM< 2000 km/s); (ii.) the forbidden line emission is relatively weak ([O III]/H/3 < 3); 
(m.) there arc often strong emission features from Fe II and high ionization optical lines. It 
has recently been discovered that they also have distinctive X-ray properties. ROSAT PSPC 
observations found that the soft X-ray spectra are systematically steeper than "classical" Seyfert 
Is and that the photon index appears correlated with the optical line width. NLSls also very 
frequently exhibit rapid and/or high amplitude X-ray variability (Boiler, Brandt & Fink 1996; 
Forster k. Halpern 1996 and references therein). 

NLSl observations with ASCA, which has better energy resolution and a larger band pass, 
find that the steep spectrum in the soft X-ray band is primarily due to a strong soft excess 
component with characteristic blackbody temperature in the range 0.1-0.2 keV and a relatively 
weak hard power law. The hard X-ray power law slope is either remarkably variable (Leighly 
et al. 1996; Guainazzi et al. 1996), or significantly steeper than found in broad-line Scyfert Is 
(Pounds, Done & Osborne 1995; Brandt, Mathur & Elvis 1997). The combination of strong soft 
excess and steep power law prompted Pounds, Done k, Osborne (1995) to postulate that NLSls 
represent the supermassive black hole analog of Galactic black hole candidates in the high state. 

2. Data and Analysis 

We considered all the ASCA data from NLSls in the archive as of February 1997, and three 
proprietary data sets, yielding a sample of 16 objects. A standard uniform analysis was applied to 
all data (e.g. Nandra et al. 1997; details in Leighly et al. 1997). 

In most cases, a power law model fit the spectra poorly, and an additional soft excess model 
was required. The soft component in these NLSls is hot and prominent, and a single black body 
is not broad enough to fit the whole energy range from 0.48 kcV. A disk blackbody model, which 
is the sum of blackbodies (Makishima et al. 1986), is broader. Therefore, we fit a disk blackbody 
model in the energy band > 0.48 keV. However, our discussion of absorption features below does 
not depend critically on the assumed model of the soft continuum. Fits over a truncated energy 
band > 0.6 keV with a single blackbody give consistent parameters for the absorption features, 
with comparable x^- 

In many cases, the power-law plus disk blackbody model gave an adequate fit. However, a 

significant dip remained near 1 keV in the fit residuals from IH 0707-495, IRAS 13224-3809, and 
PG 1404-F226 (hereafter H0707, IR1322, and PG1404; Fig. 1; also see Hayashida 1996; Otani, Kii 
&L Miya 1996; Comastri, Molendi & Ulrich 1997). Similar features have been reported in ROSAT 
spectra from Akn 564 (Brandt et al. 1994) and PG1404 (Ulrich-Demoulin & Molendi 1996). 
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Fig. 1. — Ratio of data to a power-law plus blackbody model, showing the edge- like feature around 
1 keV in these three objects. Filled and open symbols denote summed SIS and GIS spectra, 
respectively. 

Addition of an edge to the model improved the fit by Ax^/degrees of freedom (d.o.f.)= 33/265, 
46/307 and 35/156, significant at > 99.9% confidence level, for H0707, IR1322 and PG1404, 

respectively. For H0707 and IR1322, an additional edge reduced the by 24 and 11, respectively. 
The best fit parameters for these two-edge fits, and the single edge fit for PG1404, are listed in 
Table 1. 

Absorption edges, a signature of highly ionized gas in the line-of-sight, are commonly found 
in the X-ray spectra from bright Seyfert 1 galaxies (e.g. Reynolds 1997). Oxygen edges from 
O VII and O VIII at 0.74 and 0.87 keV respectively are expected to dominate the absorption 
profile because of the large abundance and cross sections. However, the edges in our objects are 
found at much higher energies than could be attributed to oxygen in the AGN rest frame. In this 
energy range, absorption by neon and iron is expected, but for abundances near solar, it is difficult 
to produce deep absorption edges from these elements without also finding strong oxygen edges. 
Since neon differs in atomic number by only two from oxygen, its ionization state should not be 
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Fig. 2. — When the edges and hnes are identified with transitions of oxygen, the ejection velocities 
from the nucleus can be found, contours show 67%, 90% and 99% confidence intervals; for 
IR1322, only 67% and 90% confidence are shown, a.) Absorption edge fits; b.) Absorption line 
fits. Solid, dashed and dotted lines indicate O VII, O VIII, and N VI, respectively. 



drastically different. Since both oxygen and neon are plausibly produced predominately by Type 
II supernovae, their relative abundance should not vary much. Iron species which could absorb 
from the L shell are still present when oxygen is nearly fully ionized, but the expected dominant 
ions have higher energy absorption edges than we find here, near 1.2-1.4 keV. 

We can attribute the absorption edges to ionized oxygen if the absorbing material is being 
accelerated away from the nucleus, similar to broad absorption line galaxies (BALQSO; see 
Weymann 1995 and Turnshek 1995 for recent reviews). If we identify the observed edges in H0707 
and IR1322 with O VII and O VIII, we infer an ejection velocity of ~ - 0.2 to -0.3c (Fig. 2a). 
This interpretation is supported by the consistency in velocity for both ions. For PG1404, we 
detect only one edge, and infer velocity of either —0.35 or —0.20 c depending on whether it is due 
to O VII or O VIII. 

Alternatively, the absorption features may be due to resonance line absorption. The ASCA 
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SIS detectors have moderate resolution 60 eV at 1 keV), and since the hnes have low equivalent 

width ~ 50 eV, they can be detected only if the velocity dispersion is (T^e/>3000 km s'^. We 
added "narrow" (a = 0) gaussian lines to the power law plus disk blackbody model, effectively 
modeling features narrower than the detector response width. Addition of an absorption line to 
the disk blackbody plus power law model reduced the by 30/265, 25/303, and 21/151 d.o.f., 
and the second line by 15/263, 6/301, 14/149 d.o.f. for H0707, IR1322 and PG1404 respectively. 
Addition of a third line for H0707 and IR1322 reduced the by 20/261 and 5/298 d.o.f. Details 
of these multiple-line fits are shown in Table 1. The final a-^e comparable to the edge fit results 
so we cannot distinguish these models statistically. Note that the two lines present in all spectra 
are consistent in energy from object to object, although the differences in cosmological redshift 
{z = 0.041, 0.067 & 0.098 for H0707, IR1322 and PG1404, respectively) imply a difference of 50 eV 
in observed energy. Assuming the absorber is completely opaque and fitting with the XSPEC 
notch model (equivalent to a very saturated absorption line) gives a velocity profile FWHM of 
~ 10 — 14, OOOkm/s, similar to that found in BALQSOs. Resonance absorption is expected to be 
strong in hydrogen-like and helium-like atoms, and the ratios of the energies should be the same 
for every element. Again, we can plausibly identify the two common lines in all spectra as O VII 
and O VIII at ~ 0.565 and 0.651 keV, respectively. A lower energy line in the H0707 spectrum 
is plausibly N VII at 0.498 keV, while a higher energy line in the IR1322 spectrum has no clear 
identification. These line identifications imply consistent ejection velocities, near —0.57 c (Fig. 2b). 
Absorption edges should accompany the resonance absorption lines (Madejski ct al. 1993). Adding 
O VII and O VIII absorption edges to the model with energies fixed at the values predicted if 
v/c = —0.57, we find that the optical depths r are consistent with zero, and the Ax^ = 4.61 upper 
limits were 0.27, 0.68 and 0.50 for O VII and 0.32, 0.13 and 0.26 for O VIII, for H0707, IR1322 
and PG1404 respectively. In Section 3 we show that estimations of the column densities derived 
from these upper limits are consistent with those from the absorption line equivalent widths. 

As noted above, the measured edge energies seem inconsistent with those expected from 
iron-L species. However, iron is a complex ion potentially producing a complex absorption profile 
which might fit these moderate resolution spectra. We used a power law plus disk blackbody 
model, transmitted through an ionized absorber in the AGN rest frame with iron abundance free 
(Magdziarz & Zdziarski 1995). In this model, the electron temperature T and the ionization 
parameter ^ are both free parameters, and we chose two representative values Tiow = 3 x 10^ K and 
Thigh = 3 X 10^ K. These models produced a worse fit overall compared with the models described 
above and required rather large iron overabundances. For T = Tiow, we obtained iron abundance 
lower limits of 15, 29 and 17 times solar with AxVd-o.f. = ll/264, 13/301, & 6/150, and for 
T = Thigh, we obtained iron abundance best fits values of 3.9, 7.3 and 12.8 with Ax^/d.o.f.=0/264, 
12/301, & 7/150 for H0707, IR1322 and PG1404 respectively. Iron enhancements up to 10 with 
respect to oxygen might be expected in some cases, but probably not as large as 20 (Hamann & 
Ferland 1993). A high iron abundance could not be responsible for the strong Fe II emission in 
NLSl, because the resulting cooler temperature decreases the number of Fe^ ions while the optical 
depth to escaping Fe II increases (e.g. Joly 1993). 
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In summary, we tried to explain the 1 keV features using three different models. Either 
absorption edges or lines produced a good fit which could not be distinguished statistically. The 
iron overabundance model generally gave a poorer fit. 

3. Discussion 

The outflow velocities of 0.2 — 0.6 c inferred in these objects are very large, larger than those 
observed the UV spectra of BALQSOs. Such large velocities could be difficult to identify in the 
UV, because the prominent C IV line would be shifted out of the band pass or into the Lya forest 
where it might be hard to distinguish. Also, very large velocity dispersions with moderate column 
densities might result in lines so broad that they blend in with the continuum. 

It is beyond the scope of this paper to speculate on the mechanism required to accelerate 
material to these very large velocities, but it is interesting to note that the velocity implied by the 
edge fits is close to that seen in the Galactic jet object SS 433 of 0.26 c and also to the terminal 
velocity predicted by "line-locking" of 0.28 c (Shapiro, Milgrom & Rees 1986). The larger velocities 
implied by the absorption line fits is intriguingly close to the energy of a particle in circular orbits 
around a Kerr black hole (e.g. Shapiro Sz. Teukolsky 1993). However, it seems difficult to relate 
this fact to a physical outflow mechanism. 

We obtain lower limits on the equivalent hydrogen column. For the absorption edge model, 
assuming cross sections of 2.8 and 0.98 x 10~^^ cm~^ for O VII and O VIII, respectively, and an 
oxygen abundance relative to hydrogen of 8.51 x 10~^, the O VII+0 VIII equivalent hydrogen 
column densities arc in the range 0.4 — 1.3 x 10^^ cm^^. For the absorption lines, assumed to be on 
the linear part of the curve of growth, oscillator strengths of 0.694 and 0.416 for O VII and O VIII 
were used, giving O VII+0 VIII equivalent hydrogen column densities of 1.6 — 2.1 x 10^^ cm~^. In 
each case, the column density upper limits on the absorption edges predicted to accompany these 
absorption lines were larger, showing that this model is viable. While estimation of the ionization 
parameter depends on the input continuum and is beyond the scope of this paper, we note that 
the O VIII column is always larger than the O VII column, implying a fairly high ionization 
parameter. 

Rest frame absorption features in the X-ray spectra of broad-line hard X-ray selected AGN 
are common (Reynolds 1997), plausibly arising in the same material as ~ Zab absorption lines 
found in the UV (e.g. Mathur, Elvis &; Wilkes 1995). These 'associated' absorption features may 
be related to the broad UV absorption lines seen in higher luminosity objects (e.g. Kolman et 
al. 1993). Evidence suggests that some aspect of the NLSl central engine is significantly different 
compared with broad-line objects; for example, they may be characterized by a higher accretion 
rate relative to Eddington (Pounds, Done &; Osborne 1995). The rapid, higher amplitude and 
perhaps characteristically nonlinear X-ray variability may be evidence for strong relativistic effects 
(Boiler et al. 1997; Leighly et al. 1997). These results may indicate a higher level of activity 
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relative to the black hole mass in NLSls, and strong relativistic outflows might be expected. 

The blue-shifted absorption features discussed here are reminiscent of those found in the UV 
spectra of broad absorption line quasars. A connection between NLSls and BALQSOs may be 
quite reasonable considering that they have many optical emission line and broad band continuum 
properties in common. Many NLSls and BALQSOs show strong or extreme Fe II emission and 
weak [O III] emission. Objects which have the strongest Fe II emission and weakest or no [O III] 
emission tend to be either low ionization BALQSOs or NLSls (Boroson &; Meyers 1992; Turnshek 
et al. 1997; Lawrence et al. 1997). Many NLSls and low ionization BALQSOs have red optical 
continuum spectra, and relatively strong infrared emission (Boroson k, Meyers 1992; Moran, 
Halpern & Helfand 1996; Turnshek 1997). Finally, both classes are predominantly radio quiet 
(Stocke et al. 1992; Ulvestad, Antonucci, & Goodrich 1995). 

An interesting possibility is that the low- ionization BALQSOs and NLSls have a common 
parent population (e.g. Lawrence et al. 1997). If so, perhaps objects with intermediate properties 
between NLSls and low-ionization BALQSOs should exist. It was recently reported that NLSl 
IRAS 13349+2438 has UV broad absorption lines (Turnshek 1997). But while most BALQSOs are 
X-ray quiet, this object is a bright soft X-ray source and has the very steep hard X-ray spectrum 
and rapid X-ray variability characteristic of NLSls (Brinkmann et al. 1996). 

KML acknowledges many enlightening discussions with Jules Halpern, and helpful advice 
from Tim Kallman. KML gratefully acknowledges support through NAG5-3307 {ASCA). KN 
thanks the NRC for support. 
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Table 1. Spectral Fitting Results 



Target 




Power-law 


Blackbody 




Absorption Feature 








Index 


Norm.'^ 


kT (keV) 


Norm."^ 


Energy (keV) 


Depth 


ID« 










ABSORPTION EDGE MODEL: 








H0707 




^•^^-0.21 




20+°-''2 
^■^"-0.03 


26+°-^'' 
^■^"-0.09 


1.09 ± 0.03 


0.57 ±0.20 
0.67 ±0.23 


VII 
VIII 


318/263 


IR1322 
PG1404 


0.48l°-«3 

U.OO_(j 34 


1 97+0-26 
-0.25 

1.67 ±0.42 




0-191^3 


r, 1O+0.49 

U.iO_o.07 


1.03 ± 0.04 
1.07 ±0.03 


U.DU_Q 25 


VII 
VIII 

VII or VIII 


339/300 
153/151 








ABSORPTION LINE MODEL: 








H0707 


1.6 ±0.08 


2.25 ±0.18 




0.15 ±0.02 




0.98 ± 0.03 
1 19+0.03 
^•-^^-0.02 
1.24 ±0.03 


47+18 
46 ±20 


N VII 
VII 
VIII 


319/261 


IR1322 


2 S+l-l 
^••^-0.9 


1.89 ± 0.23 




0.13 ±0.02 


t: 0+22.2 
0-i_4.1 


1 09+0-05 
1.23 ± 0.04 


24+18 

^^-17 

^°-24 


VII 
VIII 


343/298 


PG1404 




1.68 ± 0.33 


1 t:o+0.75 
^■^•^-0.54 


1O+0.03 


4 8+34 


1 00+0. 06 
1.00_Q 05 

1.12 ±0.03 
1.24 ± 0.04 


461-^ 
5111? 


? 

VII 
VIII 


150/149 



Note. — The quoted uncertainties arc 90% confidence for two parameters of interest (Ax^=4.61). " Neutral continuum aborption, in 
units of 10^1 cm^2. The lower limit is the Galactic absorption column for the edge fits. Absorption feature depths arc the optical depth 
r and the line equivalent width in eV for the edge and line models, respectively. Listed on successive lines are parameters for models with 
more than 1 absorption feature. ° In units of 10"* photons cm"^ s~i at 1 keV. Multiple black body model (MaJjishima et al. 1986) 
xlO^ {Ri N / Dio)^ * cos{6), where Rin is the inner disk radius in kilometers, D is the distance to the source in units of 10 kiloparsecs, and 
is the disk inclination. " Tentative edge/line identification, implying velocities in Fig. 2. 



